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^ . Abstract 

m 
o 

o . 

' totally asymmetric simple exclusion process (TASEP) with a blockage that lowers the jump rate 

CO ' at one site from one to r < 1. Exact finite volume results serve to bound the allowed values for the 



We present new results for the current as a function of transmission rate in the one dimensional 



current in the infinite system. This proves the existence of a gap in allowed density corresponding 
to a nonequilibrium "phase transition" in the infinite system. A series expansion in r, derived from 
the finite systems, is proven to be asymptotic for all sufficiently large systems. Pade approximants 
^ ' based on this series, which make specific assumptions about the nature of the singularity at r = f , 

^ I match numerical data for the "infinite" system to a part in f O**. 



^ ! 1 Introduction 

The one dimensional totally asymmetric simple exclusion process (TASEP) is a continuous-time sto- 
chastic process in which particles on a one dimensional lattice jump independently and randomly 
at unit rate to vacant neighboring sites on their immediate right It corresponds to a Kawasaki 
exchange dynamics at infinite temperature and infinite electric field The stationary state of 
this system for N particles on a ring of K sites, K > N, gives equal weight to all (^) permissible 
configurations. This measure goes over, in the limit K — > 00, N/K p, to the product measure 
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with occupation probability p. The TASEP is thus the simplest driven diffusive lattice-gas model 
whose dynamics does not satisfy detailed balance Q. It is also, for the infinite lattice, an example 
of a microscopic system from which one can derive Euler-like hydrodynamical equations [^, e.g. the 
Burgers equation. 

In an earlier work we introduced a variant of the TASEP where the jump rate across one bond 
of the system was reduced from ltor,0<r<l. If one thinks of the TASEP as a model for fluid flow 
in a pipe, this is analogous to a restriction in the diameter of the pipe. More realistically perhaps we 
could consider a superionic conductor @ like Agl or KAg4l5, in the geometry of a pinched doughnut, 
with a time-varying magnetic field generating an electromotive force; or a road under construction in 
a model of traffic flow . Clearly the rate decrease will increase the particle density to the immediate 
left of this "blockage" bond and decrease the density to its immediate right, but what is not obvious is 
that this perturbation may have global in addition to local effects. Equivalent to introducing a defect 
into a growth surface this blockage can cause the nonequilibrium stationary states of the model 

to exhibit a segregation into high- and low-density regions. This allows the full complexity of the 
model, previously available only through time-dependent studies, to be displayed via the stationary 
state |5|, |lO|. Among other features we observed for the system of N particles on a ring of size K, 
with N and K large, was the equivalence of the exponent governing the time and space scaling of 
shock fluctuations. 

It is convenient to label the sites on the ring from —K/2 to K/2 {—K/2 and K/2 refer to the same 
site), with the blockage located at the bond between and 1. The stationary density profile for a 
half-filled system can be seen in fig. As we approach the limit K ^ co keeping N/K = p fixed, we 
have that far from the blockage the density becomes uniform and the stationary state appears to be 
asymptotically a product measure [^] with densities lim^_^-|-oo p{x) = p±, with pj^ < p^, independent 
of /9, for 1/9 — 1/2 1 less than a certain value described in the next section. The densities are related by 
the condition of constant current: J = — p+) = — /)_). When segregation takes place we 
thus have 1/2 — = — 1/2. One is of course interested in the asymptotic densities p± and current 
J as a function of the transmission rate r. In these could only be determined numerically; the 
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Figure 1: Density profile (from time average in simulation) for a half- filled system with 600 sites and 
periodic boundary conditions. The blockage bond is located between site and site 1 and has the 
value r = 0.33. 



simple estimate obtained by neglecting the (large) correlations near the blockage, i.e. by assuming a 
uniform product measure with density p- (resp. p+) to the left (resp. right) of the origin, gives 

J rp_(l -/)+)« (1) 
(1 + rj^ 

which is accurate only for very small r. This is in contrast with the translation invariant case, where, 
as already mentioned, there are no correlations in the stationary state. Here, in order to obtain more 
accurate estimates, it is necessary to include the effects of correlations, which decay slowly with the 
distance from the blockage, typical of systems with conservative dynamics not satisfying detailed 
balance |l|, |l|, 0]. (In § this decay was found numerically to go like the inverse power of the 
distance.) 

Although one is interested in asymptotic behavior, it can be fruitful to examine the behavior 
of small systems with open boundary conditions for which one can determine the stationary state 
exactly. In a closely related model this led to an exact solution for all system sizes. Furthermore, 
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as we shall see, results for such systems give a systematic improvement on the estimate (|Tj) as well 
as bounds on the current J for the infinite system at fixed r and thus also on the maximum value 
of p+. We therefore studied systems, ranging in size from 2 to 10 sites, in which particles are added 
at the left with rate a and removed at the right with rate /?. At non-boundary sites particles jump 
independently to empty neighboring sites on the right with rate r for the jump between sites and 
1, where the blockage is located, and with rate 1 for all other jumps. The process is thus defined 
(on a lattice with K = 2L sites) by the generator C giving the rate of change of any function of the 
configuration rj = {r](—L + l),r]{—L + 2), . . . ,ry(L)}, where i]{k) = or 1 is the occupation number 
at site k: 

-1 

^fiv) = a[/(r/_L+i,+) - /(r/)] + E [fiv^^^+l)-f{v)] 

i=-L+l 
L-1 

+r[f iVo^i) - fiv)] + E - fiv)] (2) 

i=l 

where rji^j{k) gives the configuration at site k after an attempted jump from site i to site j: 



r]i^j{k) = < 



the boundary terms are given in terms of 



1, k = j and r/(i) = 1, 

0, k = i and r]{j) = 0, (3) 
r/(/c), otherwise; 



m,±{k) = { (4) 

rj{k), k ^ i. 

2 Exact Solutions For Small Systems 

For a system with 2L sites there are 2^^ possible configurations. By considering the set of equations 
{Cfi) = where fi is the characteristic function of the ith configuration, one obtains a system of 
equations whose (normalized) solution is the unique stationary state for the model, i.e. one obtains 
the unique stationary solution of the master equation for the open finite system. 



Solving such a system of equations is a daunting computational task. However, by fixing a and /3 
and then making use of the reflection symmetry of the problem (which for a = /3 reduces the number 
of variables by a factor slightly less than 2) we were able to exactly solve for systems up to size 8 
leaving r as an indeterminate parameter and size 10 for fixed r. The computations were performed 
using Maple V running on a SPARCstation 2. 

Typically we took a = [5 = 1; the asymptotic behavior should be independent of a and (3 provided 
they are greater than some critical value which goes to zero as r goes to zero and equals 1/2 for r > 1 . 
This critical value is simply related to the maximum current the infinite system with blockage r can 
support; it is just the density of a product measure at the maximum current: (1 — \/l — 4Jmax) /2. 
Clearly Jmax(?') is a monotone nondecreasing function of r equal to 1/4 for r > 1. The system 



"selects" the state of maximum current, p+ = (1 — \J\ — 4Jinax)/2 and p_ = 1 — if the boundaries 
can supply and remove particles quickly enough [|l^] , so the (asymptotic) state should not depend on 
the precise values of a and /?; corrections to this asymptotic behavior will be local to the boundary 



region [15, 17|. Alternatively, if one considered periodic boundary conditions instead of an open 
system we would expect equivalent asymptotics with densities /9_ = 1 — />+ for all values of the 
average density p between (1 — \/l — 4Jinax)/2 and (1 + ^Jl — 4Jmax) /2. For average densities outside 
this interval the limiting asymptotic densities should be equal to the average density both to the left 
and right of the blockage (see fig. ^) . 

We present some of the results for small systems below. JL(r, a) is the steady state current in a 
system with 2L sites at a = (3 = a; Jiir) = JL(r, 1). We have 

J^{r,a) = , Ji(r) = , Ji(r,l/2) = ; (5) 

^ ' 2a + 3r' ^ ' 2 + 3r' , / ; 1 + 3^' ^ ' 



J2 {r, a) = 2ar 



Aa^ + 8a^ + 4o^ + (ja + Ua'^ + Ga^) r + (2 + 4a + 20^) r^ 

8a^ + 16a^ + 8a^ + (Ua^ + 36a^ + 36a^ + 12a^) r (6) 

+ (l8a + 45a2 + 440^ + 18a^) + (5 + 14a + 14a^ + 6a^) , 

J (^) _ + 25r + 8r^) 2r(9 + 29r + 18r^) 

^^^^ ~ 32 + 98r + 125r2 + 39r3' ^ ~ 18 + 85r + 215r2 + 130r3 ' ^' 



J3(r) = 4r (47775744 + 261095424r + 669424384r^ + 1009680576r^ + 968982368r^ 
+609395274r^ + 250834237r^ + 65287925r^ + 9784215r^ + 644781r^) / 
(191102976 + 1331036160r + 4447316992r2 + 9277942272r3 

+12731145304r^ + 11671707972r^ + 7170513506r^ + 2914237861r^ 
+753023405r'^ + 112354075r^ + 7383541ri°). 



(8) 



We observe that Jl{t) is a rational function of r with integer coefficients, and the order of the 
function and the complexity of the coefficients grow rapidly with system size — J^ir) has terms up to 



34th order with the largest coefficient being 51 digits long, and is reproduced in Appendix A. We 



note that the value of Jl(1) was obtained explicitly in |15| for all L, namely that 



-^^(1) = 7 + 



4 4(1 + 4L) 



Joo(l) 



1 + 



1 + 4L 



(9) 



3 Bounds On The Infinite System 

For r = 1, product measure with any density p £ [0, 1] is stationary in the infinite system, and thus the 
system can have any current in the range to 1/4. For r < 1, this is not necessarily the case; Jmax(^) 
may be less than 1/4. In this case only a range of densities satisfying \p± — 1/2| > \/l — 4Jinax(?')/2 
is permitted. The problem thus is to find bounds on Jmax{r). A very simple bound on Jmax(?') can be 
obtained by noting that if we remove the right (or left) half of the system we are left with a system 
of L sites with input (removal) rate a and removal (input) rate r. Calling J{L; a, (3) the current in a 
system of L sites with input rate a, removal rate /3, and all "internal" jump rates now being unity, 
we clearly have 

Jdr, a) < J{L; a, r) = J{L; r, a), (10) 

where the last equality, as well as exact formulae for J(L;a,(3), have been computed in |jT5|, |l^. In 
particular for a > r and r < 1/2 we have Uuil^oo JiL; a, r) = r(l — r) so that Jmax('^) ^ t{1 — r) for 
r < 1/2. 
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Figure 2: Jiir, 1) and Jl(?', 1/2) for several values of L. Where Jiir, 1) < 1/4 there must be a gap 
in the allowed stationary measures of the infinite system. 

For a and r both greater than 1/2 the right hand side of ([To| ) approaches 1/4 so it yields no new 
information — ^just a proof that Jmaxif) cannot exceed 1/4. 

To obtain better bounds we note that for any configuration, since the maximum rate at which a 
particle attempts to jump on to any site is bounded by one, the current in a system of size L with 
boundary conditions a = /3 = 1 cannot increase as L increases. Note that a = f3 = 1 corresponds 
to keeping the site —L always occupied and the site L + 1 always vacant in a system of size L' > L. 
Thus JLif) is monotonically decreasing in L, and for every finite value of L, Jiif) is an upper bound 
for Joo{r)- (Similar arguments show that JL{r,a) is nondecreasing in r and a.) These bounds prove 
the existence of a gap in the set of stationary measures of the infinite system; the results for finite L 
(illustrated in fig. |2|) show that wherever there exists an L such that Jiif) < 1/4 there is a forbidden 
range of currents. (One result from our simulations is that Ji5ooo(0.8) = 0.24979(5), indicating that 
the gap exists at least up to r = 0.8.) 
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We note here also that letting a — > cxd in JL{r,a) corresponds to reducing the size of the system 
from 2L to 2{L — 1) so that lima^oo Jl{t,cl) = JL-i{r). Using monotonicity in a this implies again 
that Jiif) is monotone nonincreasing in L. 

For the case a = /3 = 1/2 the sequence Jl{t, 1/2) appears to be increasing in L, and we believe 
but cannot prove that Joo{f) > 1/2)- Proving such an inequality is nontrivial since Joo(^) is the 
maximum current the infinite system can have; there exist stationary measures where the current is 
less than Jl{i~, 1/2). Thus our results for finite L only provide rigorous one-sided bounds. Note also 
that Jl(1, 1/2) = 1/4 for all L, c.f. |l|. 

Series expansion 

While there is no apparent pattern to the "raw" expressions for the current ( |A.l| ), one does 

emerge if we examine a Taylor expansion around r = 0. We obtain: 



T f ^ 3 2 9 3 27 4 81 5 243 6 7x 
JJr) = r r"^ + -r^ H r^ + 0(r') 

^ ^ 2 4 8 16 32 ^ ^ 



(11) 



So 19 , 257 4 24105 . 829297 « 7, 

-r^ + —r^ + r^ + 0(r^), (12 

2 Ifi 25fi 4n9fi 65536 ' v ;> \ j 



Joir) = r r H r r i /■ 

^ ^ 2 16 256 4096 65536 

, 3 2 19 , 21535 4 919407829 . 7398899579671 . 7, 

Jofr) = r H r^H r^ + 0(r^), (13) 

^ ^ 2 16 27648 214990848 417942208512 ^ " ^ ^ 

, , , 3 2 19 . 21535 4 77729356627 . 

J4(r) = r r"' H r"* H 

^ ^ 2 16 27648 146767085568 

1067903077191004635349 ^ ^,7, 
126214320739011526656 ^ ^ 

There is a clear pattern: as we increase the size of the system the low-order coefficients stop changing 
after a certain point. Assuming the continuation of this behavior, and including the results for the 
size 10 system, we have 



(14) 

106790.30771910046.35.349 „ 

+ - 



T . ^ 3 2 19 3 21535 4 77729356627 . , « 7, 

J^(r) = r H r'^ -\ - .3278724755(l)r^ + 0(r ) 

^ ' 2 16 27648 146767085568 ^ ' ^ ' 



(15) 



3 2 19 3 5 -59 -73 4 13 • 33613 • 177883 5 , « ^, 7, 

" 2 + 24"^ " 210-33 + 226 7^7 ^ - •3278724755(l)r6 + 0{r')- 
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in the second part of (15) we show the prime factorization of the coefficients; the denominators appear 
deceptively simple while examination of the numerators proves less instructive. 

Including the dependence of the boundary terms (i.e. taking a ^ \) does not significantly alter 
this behavior: 

Ji(r,a) = --^-' + 4^-'-^-' + 0(r^)' (16) 

, , ^ 3 2 9a^ + 18a^ + Ta^ - 8a - 7 o 
•^^(^'"^ = + 4a^(a + l)^ ' 

54a^ + 2160*^ + 300a^ + 64a'^ - 294a3 - 3700^ - 188a - 39 , ^, 

r^ + 0{r''), (17) 



16a3(a + 1)4 



3 o 19 



2 lb 



514a^'4 + 6939a^^ + 415510^^ + 144387a^^ + 316671a^° 
+ 432661a^ + 285181a^ - 176743a'^ - 702157a^ - 944908a^ - 799104a^ - 457504a^ 



172480a^ - 38528a - 3840 rV 256a-^(l + 2a)(a + 2)^(a + 1)^ +0(r^). (18) 



The dependence on the boundary appears one term earlier than in (|Tl|)-(|T4|), but otherwise the 
structure is the same. 

In fact, the assumption regarding the behavior of the Taylor coefficients for progressively larger 
systems expressed in (pT|)-(p^ can be proven to be correct: 

Theorem 1 Fix L2. Then for Li < L2, Jh^ir) = JlA^) + 0{r^^+^). 



Theorem 2 Fix L2. Then for Li < L2 and a > 0, JL2{r,a) = JlA'^) + 0{r^^+^). 

Thus we see a small system not only bounds but also provides a good approximation (at least for 
small r) for any (finite) system that is larger. This also strongly suggests that the approach of Jiif) 
to Joo{t~) is exponential for r < 1. Of course the rate of the exponential approach vanishes at r = 1, 
where the convergence becomes algebraic (see (^). We reserve the proof for Appendix B. 
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4 Pade Approximants 



The Taylor series given in (15) gives an accurate measure of the current for smah r, but for large 
r it is less successful. In fact, we can see that the series must break down by r = 1: although one 
commonly thinks of r as being a transmission rate < 1, it is perfectly acceptable to take r > 1 in the 
generator (^). It is also fairly easy to see, by comparison with a system with a boundary at the origin, 
that the infinite volume current Jooif) will have the same value for all r > 1, namely Joo(l) = 1/4, so 
that there must be a nonanalyticity for some r < 1 in Joo{f)- One can also examine the coefficients 
of ( p!5|) and see that they apparently decrease (in magnitude) rather slowly, indicating that the radius 
of convergence of the series is most probably 1. Some numerical analysis indicates that there is no 
discontinuity in any of the derivatives of Joo{f) at r = 1; this evidence leads us to hypothesize that 
there is an essential singularity at r = 1 in Joo{r)- An alternative which cannot be ruled out by our 
results is that Joo{f) = 1/4 for r > ri with ri < 1. Numerical results for the structure of Jl{t) would 
seem however to argue in favor of a changeover at r = 1. 

We thus look for a function of the appropriate form for the current. The "simplest" function with 
an essential singularity at r = 1 that also gives Joo(O) = and Joo(l) = 1/4 is a function of the form 



Joo(r) = l/4-exp[/(r)]/4, 



(19) 



where /(r) has a simple pole at r = 1. We therefore examined functions of the form 



f{r) 



1 1 



■ exp 



p{7 



[I - r)q{r) 



(20) 



specifically, we used our Taylor series to fit the following Pade approximants: 



1 1 

- - - exp 
4 4 



1 1 

4-4^"P 



1 1 
- - - exp 
4 4 



— 4r(l + air + a2r^) 
(1 - r)(l + hir + b2r'^) 

-4r(l + a'^r + a^r"^ + a'^r^) 
{1 - r){l + h'^r + b'^r'^) 

-4r(l + a'/r + ^V^) 



(1 - r)(l + 6'/r + + h'y) 



(21) 
(22) 
(23) 
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Figure 3: Taylor, Pade and numerical results. Estimated errors for the numerical calculations are 
less than the thickness of the lines, ranging from 5 X 10"^ to 2 x 10"^. 

The actual values of the coefficients of the Pade approximants are not particularly revealing, as 
one might expect. The Pade functions do, however, appear to be converging pointwise: for r G [0, 1] 
the maximum difference between any of J^^, J^^ and J^'^ is less than 2 x 10^^. 

Of course we are not interested in how well the different approximants approximate each other, 
but how well they approximate Joo{r)- We thus must compare the approximants with numerical 
simulations. We see in fig. ^ that the approximation is within the error bounds of the simulations; 
we happened to have plotted J'^^ but any of the approximants would have fit the data as well. In 
comparison, a Pade approximant that behaves quadratically or quartically at r = 1 (as opposed to 
exponentially) looks qualitatively similar but does not fall within the error bounds of the simulations. 

The simulation results were obtained by direct simulation of the TASEP dynamics; progressively 
larger systems were used until the current reached an asymptotic value. For small values of r this 
occurred quite quickly, but the needed system size grows quite rapidly as one approaches r = 1: for 
r = 0.7 we needed to investigate systems with 6400 sites and for r = 0.8 we needed to investigate 
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Figure 4: Phase diagram of allowed asymptotic density vs. transmission. The heavy line indicates 
the region excluded by currently available exact results; the thin line and the shaded region are 
determined by the Pade results. 

systems with 15000 sites. 

Our results can also be viewed as the determination of the phase diagram for the system. As 
mentioned earlier the current Jooir) is the maximum stationary current permitted by the infinite 
system — with some boundary conditions the current may be less. Accepting |T^] that the stationary 
measure is asymptotically a product measure then J^o = /Ooo(l — Poo), and the bound on the current 
is equivalent to a bound on the allowable range of densities. We plot the boundary of this range, 
the critical density, against the transmission rate in fig. ^. If the overall density of the finite system 
in a periodic box, p, was in the disallowed region, the system would have to segregate with high 
density yO_ to the left of the blockage and low density to the right of it, with /?_ + = 1 and 
p = cp^ + (1 — for some fraction c, < c < 1. 

The phase diagram is apparently quite different from that normally observed in systems that 
phase segregate, with a cusp at the critical point if one believes the Pade approximants accurately 
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describe the behavior of the system. Of course one should remember the rather unique nature of this 
type of transition. 

5 Discussion 

When a system has dynamics that do not satisfy detailed balance, the steady states do not (in general) 
have the form of Gibbs states. Thus, even a qualitative description of the nature of the steady states, 
particularly with respect to the dependence on the parameters entering the dynamics, is lacking. 
Standard perturbative techniques are generally inapplicable, and most of our knowledge comes from 
computer simulations, approximate calculations using renormalization group (universality) ideas, and 
a few exactly soluble models. 

It thus is quite remarkable that information determined from a system with 10 sites can be used to 
predict the behavior of a system with 10^ sites to within a part in 10^ for all values of the parameter 
r. This is particularly interesting given that quantities besides the asymptotic current, e.g. the local 
density near the blockage, are not given accurately by the results from the small systems. 

We hope that it will be possible to determine a general formula for the terms in our Taylor series — 
certainly the sequence of denominators (1, 2, 2*, 2^°3^, 2^^3'^) seems tractable. Since the coefficients 
are independent of system size (if the system is big enough) only reasonably-sized systems need 
to be studied in order to obtain useful results. Unfortunately direct computation of higher-order 

terms seems unlikely, given that the computational complexity grows exponentially, without further 
theoretical input. 

Even if our computations cannot be extended, there remains the possibility of proving that some of 
our qualitative description of the phase diagram is correct. Our series is only proven to be asymptotic 
for small r, but we have already seen that finite systems can serve to bound certain properties of 
the infinite system for arbitrary values of r: for example we believe that JL(r, 1/2) < Joo(r) and we 
know that Joo(r) < J^r, 1) for aU L > and all < r < 1. 
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Appendix A The function J4(r) 



Ur) = 

2r ■ (46594570041553208865114157983928672r34+ 
3272256807527488050658050196010194008r33 + 111157036726242688664509666570575952692r32 + 
2432978692972890221457116915305608070176r3i + 38563042893543905939312199308038000267376r30 + 
471629643887135562049943024025386831761624r29 + 4630830690747494363914241555175156474737952r28 + 
37501116089445029526599189649856279114367080r2'^ + 255349882255809833270215616484973179447971604r26 + 
1483034668343761760933311314169905866593307200r25 + 7427024657245470907468866632968858025038621224r24+ 
32341852979577601689520720487367994974417590384r23 + 123258959443962771434544300470844688653283339220r22 + 
413187952609538688738716350764927640900234526696r2i + 1222925851525408590841791369809880777445578325556r20+ 
3204670630666608783584495910901402180619797197008ri3 + 7449424997833123151978707539741892038741323482132ri** + 
15378015964501993765574252081802020182175729230264ri'' + 28201950376870887559972172287807374530071235163796ri^ + 
45931932106228251497346831830246712730250907133712ri5 + 66367346077080103767105589115636872559989214550441ri4+ 
84922961375392487277902239726735636329302249943128ri3 + 95985607018161363007960499193162341535912061049744ri2+ 
95495899491963130911426581603975014896548923261952r" + 83250539026363674821927230131182759997615987109888ri0+ 
63222693774680812665162035869098746973122020442112r3 + 41512233585213949665760166256158119415119899262976r8 + 
23337512260240138424809534978398421177693582131200r'' + 11088988610809857322221895274024151043294063755264r6 + 
4375523126050494396079936095098710196366170128384r5 + 1398264149827785557648870482515232428810049683456r4+ 
348505991819536772424594468421800136094590697472r3 + 63703892226452562097992685657738901012902576128r2 + 
7610495543799775598845584755474871742038540288r +446443113770610009364511999835819193606864896) 

(289620280319036422192318811402048808r35 + 20347660058076078580887803589027266100r34 + 
691509180090958636822780801273385528578r33 + 15143171306946285250897143672457580856112r32 + 
240156757626348422508288967223318773219730r3i + 2939013682631225236236144632768990651743996r30+ 
28878453872089870194234547560112273497410248r29 + 234055024231518053113128876309994685465865998r28 + 
1595219951014207692709829729475814422350403348r27 + 9274911880978847804107170124884025274187759308r26 + 
46507371429345382752005370961862741048375464600r25 + 202819905961978964692531821554688508829966150639r24+ 
774305913910573093515211395294773846677490470023r23 + 2600902978686057892600556891176773816130767065580r22 + 
7716564464804046735459659032377460737774938262294r2i + 20279584392798156542425998004346657157480554084969r20 + 
47304820202361946730446972991019456869802234988211ri3 + 98064294514560130382369885216651184747849062158254ri**+ 
180769182668553783438231100925433451981380541595764ri'' + 296287842239143809420984903211913226411657347513603ri'' + 
431494541664969696255270976652324836280276021938501ri'5 + 557613386760147331351017000992052900309321526028898ri4 + 
638162101110332154357279278430725658354907392715680ri3 + 645076256680960920649709526072469900202092833324704ri2 + 
573956167631320085801939332326120871725603509536768r" + 447552460656613237448046131746142189157548579971072ri0 + 
304176946559301647309430508519196663546007393796096r9 + 178937028617243446184563625330908654841920346914816r** + 
90290026293392171865547954006547530297965379846144r'^ + 38609332952461801047936627883941955240573198663680r6 + 
13757984465655294150768749140417907410320746348544r5 + 3986654552595495753325700874 754842100218571259904r4+ 
904824015385555468140765882234506884944960159744r3 + 151187962701066997262421713581553533038335361024r2 + 

16560320428911381225784705510457201064897675264r + 892886227541220018729023999671638387213729792) 

(A.l) 



Appendix B Proof of Theorem |T| 



To prove theorem [T|, we will need to consider systems with unequal numbers of sites to the left and to 
the right of the blockage, say L_ and L+, respectively; Jl_,l+(?^) will represent the current in such 
a system with a = P = 1. We will prove the following lemma: 



Lemma 1 For L'_ < L_ and L'^ < L+, 



JL_,L+(r) = Jl'_,l4^) + Oir'^''^^) and Jl,,lM = JL.,L'{r) + 0{r'^'++^). 
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Theorem || is a direct consequence: for L2 > Li, Jl2(^) = Jl2,L2{'^) = '^L2,ii(^) + 0{r^^~^'^) = 

Now we prove the lemma, considering only the case L'_^ < since the case L'_ < L_ is the same 
by symmetry. We will need to consider probabilities of certain collections of configurations. We write 

Pr(r/_L_+i,r/-L_+2,- • ■ ,Vo;VuV2, ■ ■ ■,Vn)l_,l'^ (B.l) 

for < < L'_^. This represents the marginal probability of the configuration at sites 
— L_ + 1, — L_ + 2, . . . A^; the rest of the sites can take arbitrary values. The semicolon indicates 
the position of the blockage in the system. 

Let eR^f]) be the number of the r]i,r]2, . . . ,1]^ that take the value 1, i.e. eR^rj) = J2iLi Vi- Fo^' 
small particles to the right of the blockage are rare and can be treated as excitations; eij(ry) can be 
thought of as the (right) excitation number in the system. 

Our proof proceeds by induction: Suppose we know all probabilities of the form 
Pr(r/_L_+i,r?_L_4.2, • • ■ ,r]N)L-,L'_^ to a certain accuracy, namely to order eR{{r]i, . . .,r]N})+ 
k in r. Then (for > 1) we can (show that we can) compute all probabilities of the form 
Pr(r/_L_+i,?7_L_+2, ... . . ■ ,r]N-i)L-,L'_^ up to order brUtji, . . . ,rjN-i}) + A; + 1 in r. 

This is sufficient to prove our lemma, and in fact more general results: for any size system (with 
L'_^ > N) we can compute Pr{r]^L-+i,r]-L-+2, ■ ■ ■ ,r]o',Vi, ■ ■ ■^Vn)l^,L'^ up to order e/j(r?) simply by 
making use of the fact that excitations are created at rate r, providing the initial step in the induction. 
This estimate (and as a result all following estimates) is independent of L'j^ and thus also valid for 
all -L_|_ > L'j^. Since the current is 

J(r) = r Pr(l; 0) = r [1 - Pr(0; 0) - Pr(l; 1) - Pr(0; 1)] (B.2) 

if we can iterate the induction step j times we know the current to order 2 + j. 

So let us return to the induction step. Consider the transitions that occur between the different 
values of r/-L_+i5 • • • Vn-i- we will give estimates on the rates of these transitions in the stationary 
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exact: 0{1) exact: 0(1) exact: 0(1) exact: 0(1) exact: 0{1) exact: 0(1) 







exact: r \ / exact: r \ / exact: r \ / exact: r exact: r 




■ ( ■ w w 




l-O(r) 1 l-O(r) 2 l-O(r) 3 1 - 0(r) 4 1 - 0(r) 1 - 0(r) A^- 1 

to order A; + 1 to order fc + 1 to order + 1 to order fc + 1 to order A: + 1 to order k + I 



Figure 5: Markov Chain-like representation for the transition rates in the partial system. The indices 
on the "states" are the eii{r]) which are altered by the "horizontal" transitions. The "circular" 
transitions are those between states with the same e/j(r/). 



state. (The rate of a transition r r', rate [r r'], where r and r' are sets of configurations, is 

rate [r r'] = ^im^(Ai)"^ ¥t{t' at time t + Atjr at time t). (B.3) 

Thus all (elementary) transitions between individual configurations occur at rate 0, 1 or r.) 
Transitions that increase eij({r/i, . . . , t^tv-i}) involve a jump across the blockage, and thus simply 
have rate r. Transitions that keep ei{({r/i, . . . , t^tv-i}) fixed do not involve sites outside of 
{— L_ + 1, — L_ + 2, . . . , — 1} and so have rates that are simple integers. Transitions that decrease 
eR({r/i, . . . , ?7jv-i}) require that we know something about site N, since the only way to reduce 
eR{{rji, . . . , ?7iv-i}) is for a particle to move from site — 1 to site N . So 



rate [ei?({r?i, . . .,m~2, 1}) ^ eR{{r]i,. . . , ?7Ar_2, 0})] 
^ Pr(77„L_+i, . . . ,ryo;m, ■ • • , W-2, 
Pr(r/_L_+i, . . . ,r/o;r?i, . . • , W-2, 

Pr(??-L_+i, ••• • -,^-2,1, 

= 1 - 



(B.4) 



Pr(?7_L_+i, . . . , r/o; . . . , r?Ar_2, 1, 1)l„,l; + Pr(r/-L_+i, ... , rjN~2, 1, 0)l_,l^ 
The numerator of (B.4) is 0{eji{{r]i, . . . ,?7Ar_i}) + l) and we know it to 0(e/j({ryi, . . . ,r/jv-i}) + A; + l) 



by the induction hypothesis. The denominator of (B.4) is 0{eji{{rii, . . . ,r/jv_i})) and we know it to 
0(e/j({?7i, . . . ,'r/Af-i}) + k) by the induction hypothesis. Thus the rate in ( [B.4[ ) is 1 — 0(r) and we 
know it to 0(e/j({?7i, . . . ,r]N-i}) + k + l). 
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Fig. ^ thus represents our transitions, where we group together states with a common cr. It is 
clear that it is consistent to solve for the probabilities in stationary state of the subsystem represented 
by fig. I to 0(eH({r/i,... , r?Af-i}) + k + 1). This completes the induction step and thus the proof of 
the lemma. 

The proof of theorem ^ proceeds in identical fashion. The only difference is that in ( |B.1| ) one is 
limited to considering < < L'^ instead of < < L'_^_, ensuring that the transition rates for 
changing the excitation number do not depend on the boundary conditions. Thus one can perform 
the induction step one time less than for the a = (3 = 1 case and the boundary dependence appears 
one term earlier in the Taylor series. 
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